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■ Abstract 

We propose a symmetry of the Dirac equation under the interchange of signs of eigenvalues 

»—5 ■ of the Dirac's K operator. We show that the only potential which obeys this requirement 

Q\ . is the Coulomb one for both vector and scalar cases. Spectrum of the Dirac Equation is ob- 

tained algebraically for arbitrary combination of Lorentz-scalar and Lorentz- vector Coulomb 
potentials using the Witten's Superalgebra approach. The results coincides with that, known 

>• ■ from the explicit solution of the Dirac equation. 

ON ' 

in : 

1 Introduction 

O; 

The Kepler problem has an additional symmetry, more precisely, additional conserved quan- 
tity, so-called Laplace-Runge-Lenz (LRL) vector (for relevant references see, for exam- 
ple, [1]). Because of this fact there appears extended algebra together with orbital mo- 
mentum, which is isomorphic to SO (4) (for negative total energies). But components of 
LRL vector as generators of this algebra do not participate in any geometric transforma- 
tions. In 1926 W.Pauli [2] considered this algebra in quantum mechanics and obtained 
Hydrogen atom spectrum by only algebraic methods. In 30-ies V.Fock [3] considered the 
^ Schroedinger equation for the Kepler problem in momentum representation and had shown 

that the Coulomb spectrum has an 5*0(4) symmetry in energy-momentum space. Afterwards 
algebraic methods attract more wide applications and it was cleared up in the last decades 
that the hidden symmetry of the Kepler problem is closely related to the supersymmetry of 
Hydrogen atom [U [5] 

The aim of our paper is to determine what the Dirac equation tells us about this problem. 
We can see that very interesting physical picture emerges. 

For the Kepler problem in the Dirac equation Johnson and Lippmann [B] published very 
brief abstract in Physical Review at 1950. They had written that there is an additional 
conserved quantity 

he 

A = a ■ fr~ l — i(—)(mc 2 )~ 1 jpi(H — mc 2 p 3 ), (1) 
e 2 

which plays the same role in Dirac equation, as the LRL vector in Schroedinger equation. 

As regards to the more detailed derivation, to the best our knowledge, is not published 
in scientific literature (one of the curious fact in the history of physics of 20th century) 
Moreover as far as commutativity of the Johnson-Lippmann(JL) operator with the Dirac 
Hamiltonian is concerned, it is usually mentioned that it can be proved by "rather tedious 
calculations." [TJ. 
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In papers [HI E], we developed rather simple and transparent way for deriving the JL 
operator. We obtained this operator and at the same time proved its commutativity with 
the Hamiltonian. 

After we consider necessity to be convinced, that the Coulomb problem is distinguishable 
in this point of view. We considered the Dirac equation in arbitrary central potential, V(r) 
and shown that the symmetry, which will be defined more precisely below, takes place only 
for Coulomb potential. 

Moreover it is remarkable that this property remains valid for arbitrary linear combination 
of Lorentz- vector and Lorentz-scalar potentials [TT]. Therefore we expect that the spectrum 
for this general case can be obtained by pure algebraic methods, i.e. without solving of 
corresponding Dirac equation. 

Our paper is organized in the following manner: First of all we define a symmetry, which 
in fact leads to Witten's superalgebra. Then we determine the general class of so-called 
K-odd operators and apply this result to pure vector and then to mixed cases. At the end 
we derive the spectrum by algebraic methods on the ground of supersymmetric approach. 
We show that the final expression for ground state energy of Hydrogen atom coincides with 
that obtained by solving of equations of motion. 

2 Survey Of Symmetry 

Let consider the Dirac Hamiltonian in arbitrary central field, V(r), which in general may be 
Lorentz-scalar or 4th component of the Lorentz vector (or their combination) 

H = a ■ p + (3m + V(r) (2) 
The so-called Dirac's ^-operator, defined as 

K = P(E ■ f+ 1) (3) 

A — > 

commutes with this Hamiltonian for arbitrary V[r), [K, H] = 0. Here £ is the spin matrix 

a0\ 




and I - orbital momentum vector-operator. 

It is evident that the eigenvalues of Hamiltonian also are labelled by eigenvalues of K . For 
example, the well-known Sommerfeld formula for the Coulombic spectrum looks like 



E n ,\k\ = m 



1/2 

(*.)' 



{n — \k\ + Jk 2 — (Za) 2 ^ 2 



It manifests explicit dependence on k , more precisely on \k\ = j+1/2. It is surprising that 
for other solvable potentials the degeneracy with respect to signs of k does not take place. 
Physically this degeneracy leads to the forbidden of the Lamb shift. Indeed, positive k = 
j + 1/2 corresponds to aligned spin j = I + 1/2, i.e. to states (si/ 2 ,P3/2, etc), while negative 
k = — (j + 1/2) corresponds to unaligned spin j = I — 1/2, i .e. to states (pi/2, ^3/2, etc). So 
the absence of the Lamb shift (si/ 2 — P1/2) is a consequence of above mentioned degeneracy 
k — > —k. Naturally for description of this degeneracy one has to find an operator that mixes 
these two signs. It is clear that such an operator, say A, must be anticommuting with K: 

{A, K} = AK + KA = (4) 
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If at the same time this operator should commute with the Hamiltonian, then it'll generate 
the symmetry of the Dirac equation. Therefore, we need an operator A with the following 
properties: 

{A, K} = and [A, H] = (5) 
It is our definition of symmetry we are looked for. 

It is worth to mention that after this operator is constructed, we will be able to define 
relativistic supercharges as follows: 

a ^ AK 
Q 1 = A, Q 2 = i-r-r 

\k\ 

It is obvious that 

{Qi,Q 2 } = 0, Q\ = Q\ 

and we can find Witten's superalgebra [15] . where Q\ = h is a so-called, Witten's Hamilto- 
nian. 



3 Symmetry Operators 

Now our goal is a construction of the operator A. We know, that there is a Dirac's 7 5 matrix, 
that anticommutes with K. What else? There is a simple theorem [9], according to 
which arbitrary (S • V) type operator, where V is a vector with respect of I and is 
perpendicular to it, (I ■ V) = (V ■ l), anticommutes with K: 

{$-V),K} = (6) 

It is evident that the class of operators anticommuting with K (so-called A'-odd operators) 
is much wider. Any operator of the form 0(£ ■ V) , where O commutes with K, but is 
otherwise arbitrary, also is a A'-odd. This fact will be used below. 



3.1 Pure Vector Potential 

Let first establish the form of symmetry operator for vector potential only V(r) = V(r). 
By using the above mentioned theorem we wish to construct the A'-odd operator A, that 
commutes with H . It is clear that there remains large freedom according to the above 
mensioned remark about O operators - one can take O into account or ignore it. 

We have the following physically interesting vectors at hand which obey the requirements 
of our theorem. They are 

r — unit radius- vector and p — linear momentum vector (7) 

Both of them are perpendicular to /. Constraints of this theorem are also satisfied by LRL 
vector A = r — ^taW xl — lxp\, but this vector is associate to the Coulomb potential. Hence 
we abstain from its consideration for now. Thus, we choose the following A'-odd terms 

(S-f), K(t-p) and Ay 5 (8) 

and let probe the combination 

A = • f) + ix 2 K(t ■ p) + ix 3 K-f 5 f(r) (9) 
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Here the coefficients Xi(i = 1, 2, 3) are chosen in such a way that A operator is Hermitian 
for arbitrary real numbers and f(r) is an arbitrary scalar function to be determined later 
from the symmetry requirements. Let's calculate 

= [A, H] — (T,- ?){x 2 V'(r) - x 3 f(r)} + (10) 

+2if3K^{— -mx 3 f(r)} 
r 

We have a diagonal matrix in the first term, while the antidiagonal matrix in the second 
one. Therefore two equations follow: 

X\ 

x 2 V'{r) = x 3 f'(r), x 3 mf(r) = — (11) 

r 

One can find from these equations: 

V(r) = ^ — (12) 
x 2 mr 

Therefore only the Coulomb potential corresponds to the above required ±k 
degeneracy. The final form of obtained A operator is the following: 

A — Xi{(E • r) — —Kit ■ p) + —K^\ (13) 
ma mr 

where unessential common factor x\ may be omitted and after using known relations for 
Dirac matrices, this expression may be reduced to the form 

A = 1 5 {a-r- — Kj 5 (H-Pm)} (14) 
ma 

Above and here a is a strengh of Coulomb potential, a = Za. Precisely this operator is 
given in Johnson's and Lippmann's abstract [6]. 

What the real physical picture is standing behind this? Taking into account the 
relation 

K(E.p) = -tp(E-±\pxT-Tx^ (15) 
one can recast our operator in the following form 

A = S • \ r — (3 \px I- fx p\ ] + —K^ 5 

\ 2ma L V mr 



One can see that in the non-relativistic limit, i.e. (5 — > 1 and 7 5 — > our operator reduces 
to 

A^A NR = a-(?-^—\pxT-Txp\\ (16) 

Note the LRL vector in the parenthesis of this equation. Therefore relativistic super- 
charge reduces to the projection of the LRL vector on the electron spin direction. 
Precisely this operator was used in the case of Pauli electron |14j . 

Because the Witten's Hamiltonian is 

and it consists only mutually commuting operators, it is possible their simultaneous diago- 
nalization and replacement by corresponding eigenvalues. For instance, the energy of ground 
state is 
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By using the ladder procedure, familiar for SUSY quantum mechanics, the Sommerfeld 
formula can be easily derived [7j. 



3.2 Inclusion of Lorentz-scalar Potential 

Let's remark that if we include the Lorentz-scalar potential as well 

H = a- p + (3m + V(r) + (3S(r) (19) 

this last Hamiltonian also commutes with .fT-operator, but does not commute with the above 
JL operator. 

On the other hand, non-relativistic quantum mechanics is indifferent with regard of the 
Lorentz variance properties of potential. Therefore it is expected that in case of scalar 
potential the description of hidden symmetry should also be possible. In other words, the 
JL operator must be generalized to this case. 

For this purpose it is necessary to increase number of i^-odd structures. One has to use 
our theorem in the part of additional O factors. 

Now let us probe the following operator 



X = Xl (E • r) + x[ [p-f)H + ix 2 K (E ■ pj + tx 3 KYf\{r) + ix' 3 Kf /3f 2 (r) (20) 

We included factor O = H in the first structure and at the same time the matrix O = (3 
in the third structure. Both of them commute with K and are permissible by our theorem. 
The form (20) is a minimal extension of the previous picture, because only the first order 
structures in f and p participate. For turning to the previous case one must take x[ = x' 3 = 
and S(r)=0. Calculation of relevant commutators gives: 

[X, H] = 7 5 (3K{^ - 2tx 3 (m + S)h(r) + ^V(r)} + (21) 
r r 

+K{t ■f){x 2 V'{r) - x 3 f[(r)} + 
+K(3 (Z-7){x 2 S'(r)-x' 3 f 2 , (r)}+ 



5 2»^(m + g) _ + s)h{r)}+ 

r 

+W J^_2z4/ 2 ( r )}(S-p) 
r 

Equating this expression to zero, we derive matrix equation, then passing to 2 x 2 representa- 
tion we must equate to zero the coefficients standing in fronts of diagonal and antidiagonal 
elements. In this way it follows equations: 

(i) From antidiagonal structures (7 5 , rf(3K) 

^-x 3 (m + 5)/i(r) + ^W(r) = (22) 

^(m + S)-(m + S)f 2 (r) = 
r 

(ii) Prom diagonal structures( K(T, ■ r), K/3{E ■ f), pK(T> ■ p) ): 

x 2 V'(r)-x 3 f[(r)=0 (23) 
°y - x 2 (m + S)V(r) + *yV(r) = 
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x 2 S'(r) - x'J' 2 (r) = 

^i-4/ 2 (r)=0 
r 

Integrating the first and third equations in (23) for vanishing boundary conditions at 
infinity, we obtain 

/i(r) = -V(r), / 2 (r) = ^S{r) (24) 
x 3 x 3 

and taking into account the last equation from (23), we have 

Mr) = 4" (25) 
x 3 r 

Therefore, according to (24), we obtain finally 

S(r) = *L (26) 
x 2 r 

So, the scalar potential must be Coulombic. 

Inserting Eq.(24) into the first equation of (22) and solving for V(r), one derives 

V(r) = — - (27) 

r x 2 (m + S) - ^ 

At last, using here the expression (26), we find 

V(r) = (28) 
x 2 mr 

Therefore we have asserted that the ±« degeneracy is a symmetry of the Dirac 
equation only for Coulomb potential (for any general combination of Lorentz 
scalar and 4th component of a Lorentz vector). 

It seems that the conservation of LRL vector is a macroscopic manifestation of symmetry, 
which is present in microworld. 

Further, if we take into account above obtained solutions, one can reduce the X operator 
to more compact form 

X = (E • r) {may + a s H) - iKj 5 (H - /3m) (29) 
where the following notations are used 

X\ x[ 
ay = , as = 

x 2 m x 2 

Here a^-s are the constants of corresponding Coulomb potentials 

rr/ \ a V Q/ x a S 

V(r) = , S(r) = 

In conclusion we want to remark, that this expression for X was obtained earlier by 
Leviatan [TT], who used the radial decomposition and separation of spherical angles in the 
Dirac equation. 

Our approach is 3-dimensional, without any referring to radial equation and, therefore is 
more systematic and rather easy. 
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4 Algebraic Derivation of Spectrum 



Now we want to obtain spectrum of the Dirac Hamiltonian pure algebraically, without any 
referring on equations of motion. Our method is based on Witten's superalgebra,established 
in Sec. II above. 

Now we explore this algebra in a manner as in paper |7j. One defines a SUSY ground 
state |0): 

h\0) = X 2 \0) = — ► X|0) = (30) 

Because X 2 is a square of Hermitian operator, it has a positive definite spectrum and one is 
competent to take zero this operator itself in ground state. By this requirement we'll obtain 
Hamiltonian in this ground state and, correspondingly, ground state energy. After that by 
well known ladder procedure we can construct the energies of all excited levels. We believe 
that this method requires more strong justification, but nevertheless we are convinced, that 
it is true. 

Let equate X = and solve H from Eq. (29): 

H = m\(a-r)a s + iK}- 1 \iK(3-a v {a-rj\ = — T N (31) 

K + a s 

where 

N =[(a ■ r)a s - iK] [iK/3 -a v (a- r)] = 

= —asav + K[Kf3 + iav{oi • r)\ — iasK(3(a ■ r) (32) 

Now we try to diagonalize this operator using Foldy-Wouthuysen [15] like transformation. 
Because the second and third terms do not commute with each others we need several (at 
least two) such transformations. 

We choose the first transformation in the following manner 

exp(zS'i) = exp ^— -/3(a ■ rjw^J (33) 



It is evident that 



Moreover 



exp(«5i)(a ■ r) exp(— iSi) = exp(2«S'i)(a • r) (34) 
exp(2Si)/3exp(— iS\) = exp(2iSi)/3 

exp(zSi)K exp(-iSi) = K, 
expitS^pKexpi-tSi) = exp^iSJpK (35) 



and 



exp(iSi)P(a ■ r) exp(— iSi) — (3(a-r) (36) 
Therefore the first transformation acts as 

N' = exp(«5i)JVexp(-i5' 1 ) = 
= —asay + K exp(2iSi) [K(3 + iay{o. ■ r)] — iasKf3(a ■ r) (37) 
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But exp(2tSi) = coshwx + t(3(a ■ r) smh.W\. Make use of this relation, we have 

exp(2iSi)[K(3 + ia v {a ■ r)] = 
f3[K cosh w\ + ay sinh wi] + K(a ■ r) [lay cosh w\ + iK sinh toj (38) 
Now in order to get rid of non-diagonal (a ■ r) terms, we must choose 

tanh^i = — — ^- (39) 
K 

Using simple trigonometric relations we arrive at 

exp(2iS 1 )[K + iay (3 ■ r)] = K~ 1 ^k 2 - a\ (40) 
Let us perform the second F.-W. transformation 



Now 



N" = exp(iS 2 )N' exp(-iS 2 ), where S 2 = -]-{a-r)w 2 (41) 

exp(iS 2 )K/3exp(-tS 2 ) = exp(2iS 2 )K/3 (42) 
exp{iS 2 )K(3{a ■ r) exp(—iS 2 ) = exp(2zS 2 )K(3(a ■ r) 
exp(2iS 2 ) = cosw 2 — i(a ■ r) sin w 2 



Therefore 



N" = -a s a v + K^k 2 - a 2 v exp(2iS 2 )(3- 
—ias exp(2iS 2 )KP(a ■ r) 



= —asay + Ky n 2 — a v [3 cosw 2 + iKy n 2 — a v P(a ■ r) sin w 2 — 

— iasK[3 (a ■ r) cos w 2 + asK[3 sin w 2 (43) 
Requiring absence of (a ■ r) terms we have 

&<? 

tanw 2 = (44) 



K 2 — ay 



Therefore 



and finally 



N" = -a s a v + K(3Jk 2 - a\ + a% (45) 



m 



H=— j{-a s a v + K(3^K 2 -a v + a 2 s } (46) 



For eigenvalues in ground state we have 



m 



-^o — — 2 
k + ai 



(47) 



—asay ± fty k 2 — ay + a 2 s 

Now let us remember the result obtained by explicit solution of the Dirac equation for 
this case [16] 
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E _ —a s a v 

m a v + (n — |ac| + 7) 5 



± 



ds a V 

\ \a v + (n - |k| + 7) 1 



+ 



(n — 


AC 


+ 7) 2 "4 


a 2 v + 


(rz 




ac 


+ 7) 2 



where 

In the ground state n = 1, j 



Er, 



m 



r 

-1/2- 



K — GU 



AC 



j + 1/2 = 1, there remains 



± 



ay + 7 2 \ \ay + 7 



+ 



7 2 ~4 
a^+7 2 



(48) 
(49) 

(50) 



which after obvious manipulations reduces to our above derived expression (47). Therefore 
by only algebraic methods we have obtained the correct expression for ground state energy. 
For obtaining of total spectrum it is sufficient now to use the Witten's algebra. Following 
the paper [7j , the ordinary ladder procedure consists in change (for our case): 



7 — > 7 + n — |ac| 

Making this, it follows from our lowest energy formula (47) the correct expression for 
total energy spectrum, Eq. (48). 



5 Conclusions 



In conclusion, by using of pure algebraic manipulations we obtained the spectrum of general- 
ized Coulomb problem of the Dirac equation for arbitrary combination of Lorentz-scalar and 
Lorentz- vector potentials. This fact demonstrates a power of symmetry considerations. It is 
interesting by itself that the requirement of validity of Witten's algebra, or equivalently, 5(2) 
supersymmetry of the Dirac Hamiltonian for arbitrary central potentials (scalar or vector) 
leads to unique role of the Coulomb potential. 
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